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Abstract 

SiC is a robust semiconductor material considered ideal for high-power application due 
to its material stability and large bulk thermal conductivity defined by the very fast 
phonons. In this paper, however, we show that both material-interface scattering and 
total-internal reflection significantly limit the SiC-nanostructure phonon transport and 
hence the heat dissipation in a typical device. For simplicity we focus on planar SiC 
nanostructures and calculate the thermal transport both parallel to the layers in a 
substrate/SiC/oxide heterostructure and across a SiC/metal gate or contact. We find 
that the phonon-interface scattering produces a heterostructure thermal conductivity 
significantly smaller than what is predicted in a traditional heat-transport calculation. 
We also document that the high-temperature heat flow across the metal/SiC interface 
is limited by total-internal reflection effects and maximizes with a small difference in 
the metal/SiC sound velocities. 
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1. Introduction 



SiC is an exciting material with possible application in high-power devices but also 
posing interesting challenges for our understanding of the necessary heat dissipation. 
There are several SiC advantages, for example, structural robustness, presence of a 
native oxide (Si02), durability, and ability to grow high-mobility SiC layers on top of 
cheap Si substrates. The set of very fast phonons, the high purity, and the long phonon 
mean free path l^fp also ensures a very large bulk thermal conductivity. On the one 
hand, this raises the promise for high-power device applications since it indicates an 
effective heat dissipation. On the other hand, the micron scale Zmfp invalidates the 
foundation of the traditional analysis based on the Fourier law of thermal transport. A 
finite-element approach that describes the local flow in terms of bulk values for the ther- 
mal conductivities breaks down because the phonon mean free path now approaches 
or exceeds typical feature sizes (thickness of conduction channel and oxide layers and 
distance between gates and contacts). We must instead turn to phonon-based calcu- 
lations to describe SiC devices that represent a natural test case for developing an 
understanding of nanostructure thermal transport P]. 

Figure 1 shows a structural schematics (upper panel) of a SiC device for which we 
document two phonon transport effects (lower panels) that significantly suppresses 
the heat dissipation. We observe that the thermal transport is completely dominated 
by the phonon dynamics even when discussing the conductance across a SiC/metal 
interface. For the phonon transmission problem we calculate the thermal conductance 
of SiC interfaces at room temperature, 6 = 300 K, and report the variation with 
general ratios of the material acoustic impedances and sound velocities. We find that 
this thermal conductance exhibits only a weak dependence on the ratio of acoustic 
impedances (which determines the long-wavelength phonon transmission) but a strong 
variation with the ratio of the material sound velocities (that define the total internal 
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reflection of phonons). In contrast to what is normally assumed, we argue that heat 
dissipation can be maximized by matching the sound velocities of the SiC and metal 
(or substrate). 

We also document a phonon Knudsen effect that causes a significant suppression of the 
in-plane thermal conductivity for nanoscale layer thicknesses at room temperatures 
and even micronscale thicknesses at liquid-nitrogen temperatures, 6* ^ 77 K. At the 
same time we show that the traditional Fourier analysis breaks down in the description 
of such small structures and stress that shifting toward a phonon based description of 
thermal dissipation will become increasingly important with a continued downscaling. 

2. Suppression of thermal transport across SiC/metal interfaces 

Figure 2 illustrates the atomistic simple-cubic (SC) model adapted to estimate the 
finite-temperature thermal transport across the SiC/metal (or SiC/Si-substrate) inter- 
face. Both the phonon frequency uj and the in-plane momentum q\\ (parallel to the 
interface) is conserved in this description of the phonon transmission. While approxi- 
mative, our model description retains the important transport consequences of phonon 
total- internal reflection j2l El • 

This SC atomistic model has previously been used to calculate the phonon transmission 
and thermal transport in Si/Ge superlattices and heterostructures [21 El- Here we 
extend the SC model to describe the binary nature of the SiC material by choosing the 
SC volume a^^ to contain exactly one physical atom, but having an average atomic mass 
Msic = (Msi + Mc)/2. The value of Osc = 2.18 A corresponds to the physical (FCC) 
lattice constant = 4.36 A of SiC |31. Within each of the materials, these SC volumes 
are linked by force constants -Fp,t;Sic with values which are fitted to the SiC sound 
velocities |11 ci^t;SiC = ctsc \/ -^p,t; SiC / ^Sic ! as expressed in the SC approximation [2J. In 
this model we also obtain acoustic impedances [21 -^/,t;SiC = Fp,t; sic^sic / c^sc- W^e use 
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a corresponding notation {Fp^t;me, ^me, • • •) to describe the SC phonon dynamics in the 
neighboring media 'me' (metal, Si, or general substrate). We assume for simplicity 
an identical lattice constant age for both the SiC and the adjacent material, and we 



approximate the elastic coupling across the interface by Kp t = Y-^p,t;SiC-^p,t;me- 

Our calculation of the phonon conductance uk across the SiC/metal or SiC/substrate 
interface follows the approach and analysis of Ref. [HI to which we refer for a more 
detailed discussion. The modes in our SC model can be polarized either parallel ('p') or 
at right angles ('t') to the interface normal. In the following we focus on the description 
of the phonon dynamics and transmission of 'p'-polarized modes. The dynamics of a 
phonon moving in the SiC or neighboring media 'me' is then unique characterized by the 



conserved set We introduce frequencies l^p, t; sic (me) = '^\/ Fp^t;SiC(me)/Msic{me) 

and the (conserved) dimensionless measure = 2 — cos(g^asc) — cos^qyttsc) of the 
in-plane dynamics. The SC-phonon dispersion relation can be expressed 



''p; SiC (me) 



1 - COs(A;siC{me)asc) 



2uj'^ — ^^t;SiC(me)0^g|| 



and serve to establish the perpendicular momentum component fcsic(me) of any phonon 
mode (tu, g||) when moving in the SiC (and/or in the metal or substrate 'me'). 

We emphasize that this SC model description includes the affects of attenuation and 
of phonon total-internal reflection [^j. The attenuation (with associated thermal- 
transport suppression) arises because there is either insufficient or excess energy avail- 
able for the perpendicular dynamics. A phonon with in-plane momentum gy > that 
approaches from the typically softer 'me' side can be blocked from entering the hard 
SiC because the mode have too little energy to sustain the in-plane dynamics. Con- 
versely, a mode (^^,^11) which propagates in the SiC (having a real value of /^sic) can 
have too much energy to enter the (softer) metal or substrate. In our SC model, both 
of these situations causes the dispersion relation (Q) to yield an imaginary value of 
fcsic(me)5 i-G-, the correct attenuated behavior on the SiC (me) side of the interface. At 



elevated temperatures and with a finite separation in the sound velocities (material 
hardness) such total-internal reflection suppresses the thermal transport contribution 
from most of the available phase space [2 EI • 

To quantify the resulting suppression of the conductance we calculate the SC estimates 
for the probability TK{uj,q\\) of phonon tunneling at general u and gy within our SC 
model. We extract Tk{uj, q\\) by solving the Newton equation of motion for the phonon 
displacement 77g||;Sic(^ < 0) and rig^^-^neis > 1) of the (atom-scale) SC elements located 
at the interface (here located between / = and s = 1). Speciflcally, we consider then 
the reflection and transmission of an incoming SiC mode 

?/g||;Sic(^ < 0) = e^^Sicasci ^ ^ g-ifcsicasJ (^3) 
//,||;me(s>l) = B e''^''-' , (3) 

and solve the harmonic equation 

/ \ ■ me r 

Vq^^;mcis = 1) [2UJ^ - nl^^ag^J /2 = ^| ; „,e (s = 2) - ?7g|| ; me(s = 1) 



+ 



Kr, 



Vgn;me{s = 1) - '7g||;Sic(^ = 0) . (4) 



The tunneling vanishes, of course, exactly when there is attenuation on either sides of 
the interface; in the case when both fcsic and /Cme are real, we determine the tunneling 
probability 

Tk{u;; ag^^) = 1 - \Af (5) 
and hence obtain a description of the transport contribution at general {Lj,q\\). 

We calculate the flnite-temperature thermal conductance across the SiC/Metal inter- 
face by use of the formal result jSJ 0] : 



ax = 

modes ' 



(6) 



de ^ 

The 'free' interface thermal conductance do that describe the phonon transmission 
across a SiC/SiC interfaces (for example, between the doped SiC substrate and the 



undoped SiC transport channel) serves as a natural reference. We typically have the 
situation axiO) ^ cro(^) at finite temperatures (for example, above the liquid- 
nitrogen temperature). 

Figure 3 shows the results for the thermal conductance, cr^, calculated at = 300 K 
for a general SiC/metal or SiC/substrate interface within the SC model description. 
The plot reports a dramatic variation of ax as a function of the ratio of sound velocities 
Cme/csic and a much weaker dependence on the ratio of acoustic impedances Z^e/Zsic- 
Inserting the value for these ratios from Table I, the figure also documents that the 
typical SiC/metal and the SiC/Si interface at = 300 K causes a dramatic suppression 
of the thermal conductance (when compared to the free conductance ctq value given 
by Cme/csic = ^me/^sic = !)• The study identify the need for further improving our 
understanding of the thermal coupling at such interfaces. 

We also stress that this suppression of uk is qualitatively and quantitatively very dif- 
ferent from the traditional conductance variation in low-temperature Kapitza effect. In 
the low-temperature regime the transmission is given by classical acoustics and hence 
exclusively defined by the ratio Zmc/Zsic- In the present finite-temperature regime 
it is instead the total-internal phonon refiection and hence Cme/csic that specifies the 
dominant thermal-conductance variation. In general the contact region between SiC 
and metals possesses a complex structure. Our calculation of the ax variation assumes 
a perfect junction between the SiC and the metal or Si substrate, rests on a simplified 
description of the characteristic phonon modes, and can thus only provide an approxi- 
mative description. Nevertheless, our results strongly indicates that it is the difference 
in sound velocities rather than (as often assumed) the ratio of acoustic impedances 
which should be matched to maximize the high-temperature thermal transport out 
from a SiC power devices at contacts, gates, and cooling fringes. 
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3. Suppression of the thermal transport in SiC layers 



We calculate the phonon Knudsen effect on the in-plane thermal transport by solving 
the linearized Boltzmann transport equation (1-BTE) in the relaxation-time approxima- 
tion [B]. We consider a layer of thickness c^iayer in which a small thermal gradient, V0||x, 
induces a change 5N from the equilibrium phonon distribution Nq = i/f^g^'^/'^s^ _ 
The upper, lower materials boundaries are placed at 2;>,< = ±(iiayer/2 and we here use 
^2/^11 to denote the component of the total phonon momentum, q, that is perpendicular 
to the layers/parallel to the in-plane temperature gradient V^, respectively^. The bulk 
phonon mean free path, /^fp = v^r^, is given by a product of the phonon group velocity 
and the (bulk) relaxation time. We assume, for simplicity, that /^fp is independent of 
the mode and momentum and we fit /mfp against measurements of the bulk-SiC thermal 
conductivity [71. Solving then the 1-BTE 



[ q dz /mfp J q\d9 
in the absence and presence of the first term, we obtain and compare the distribution 
changes, 5A^buik and 

(^iVlayer = ^^VBulk (1 - hq^{z, p+, . (8) 

These distribution changes, in turn, determine the phonon thermal conductivity 



«:bulk,layer = -(V^)'^ ^ / ^ / (^^^.^^ 5iVb.ik,layer (9) 

in the bulk and in the layered structure, respectively. 
The interface scattering causes a suppression 

h,^{z, p+, p.) = Y: p_)e-l^-^^l/{'-p(l«-l/'?)> (10) 

which depends on the fraction of specular and diffusive 1 — scattering at the 

upper (lower) material boundary. We determine the amplitudes 0^^'^ of this variation 

^The other in-plane momentum component is denoted qy\ we seek to keep a consistent notation 
between the two calculations although the direction of the temperature gradient differs. 



7 



from an analysis of the phonon scattering at the two interfaces, z>^< = ±(iiayer/2, 
following Refs. [HIHlin]- Specifically, for phonons that travel toward = d/2 (with 
perpendicular momentum > 0), the appropriate boundary condition for specular 
scattering is 6N{Q, q^) = 5N{VL, —qz), while the corresponding boundary condition for 
diffusive scattering is 6N{Q, < 0) = 0. A similar treatment yields the boundary 
conditions at the lower interface, z<, and an approach to determine the full variation 
of the suppression (fTUj) by interface scattering. 

We introduce p = (p+ + p_)/2 and p = y/p^p_ and express the resulting phonon 
Knudsen effect 

^iv, p^, p_) = 1 _ A I (1 _ p) _ 4 f Q(j, r^, p, p) I (11) 
using special functions 



V, P, P) = iP' - 2p + l)f-'Dijr]) + {p-p){p + 1) V(^-i)D(2jr7) (12) 

and 



The analytical result extends earlier results [71 E] and permits an efficient evaluation 
in the case where the degree of specular and diffusive scattering differs at the upper 
and lower boundary. The result (jl3|) describes, in particular, the case when a thin 
(undoped) SiC transport channel is located below an oxide (p+ = 0) and on top of a 
doped SiC substrate (p_ ^ 1) or on a Si substrate (0 < p_ < 1, depending on the 
quality of the interface). 

Figure 4 documents the significant phonon Knudsen effect caused by the interface 
scattering in two typical structures (inserts) at two different temperatures. The figure 
shows the variation in the effective thermal conductivity Keg evaluated as the weighted 
average of the (interface-limited) thermal conductivities describing the (finite) oxide, 
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(finite) SiC and (finite) Si-substrate layers. The figure documents that the Knudsen ef- 
fects causes a very significant suppression of the thermal transport even up to hundreds 
of nanometers (micrometers) in the room (liquid-nitrogen) temperature cases. 

Figure 4 further documents that calculations based on the traditional approach using 
Fourier law of heat conductions is inadequate due to the micron-scale phonon mean 
free path. The panels compares calculations with such a finite-element type description 
(dotted curves) against the present phonon-based transport calculations assuming both 
diffusive (sohd curves) and specular (dashed curves) scattering at the SiC/Si-substrate 
interface. We find that the traditional approach completely breaks down at lower 
temperatures, and that moving to phonon-based calculations of thermal transport is 
essential for quantitative discussions of heat dissipation at room temperatures. 

4. Conclusions 

In this paper we have documented significant suppressions of the thermal transport 
both across a SiC/metal or SiC/substrate interfaces and within the SiC layer between 
a top oxide and the substrate. For the room-temperature thermal conductance we find 
that it is the ratio of material sound velocities, rather than of the acoustic impedances, 
which should be adjusted to maximize the thermal transport. For the thermal transport 
in a layered structure we further documented that interface scattering causes a signifi- 
cant suppression for nanoscale (micronscale) layer thicknesses at room (liquid-nitrogen) 
temperatures. Neither of these heat-transport effects can be addressed within the tra- 
ditional description of heat transport based on the Fourier law. We conclude that the 
introduction of multiple interfaces and/or continued reduction of the thickness of the 
substrate or the conduction channels has adverse effects on the overall device heat dis- 
sipation. The ongoing quest for miniaturization motivates a shift toward phonon-based 
calculations of the thermal transport in nanostructure devices. 
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TABLES AND TABLE CAPTIONS 



Table 1: Materials properties and model parameters for the calculation of the 
room-temperature thermal transport across SiC/metal and/or SiC/Si interfaces. 



«o M a.e Q Q (fe), (fe)^ (^)^ (^)^ 

A a.m.u A (m/s) (m/s) 



SiC(FCC) 


4.36 


20.05 


2.188 


13035 


6257 


1 


1 


1 


1 


Si(FCC) 


5.43 


28.09 


2.715 


8470 


5340 


0.65 


0.85 


0.48 


0.63 


Al(FCC) 


4.05 


26.99 


2.551 


6360 


3130 


0.41 


0.42 


0.49 


0.50 


Ni(FCC) 


3.52 


58.71 


2.218 


5810 


3080 


1.25 


1.38 


0.45 


0.49 


Cu(FCC) 


3.61 


63.55 


2.274 


4760 


2320 


1.02 


1.04 


0.36 


0.37 


Ag(FCC) 


4.09 


107.87 


2.576 


3640 


1690 


0.92 


0.89 


0.28 


0.27 
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FIGURES AND FIGURE CAPTIONS 
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Figure 1: Structural schematics (top panel) and summary of key phonon transport 
effects (bottom panels) in a typical nanoscale SiC device. Oxide and substrate layers 
together with SiC/metal interfaces surrounds the thin (undoped) SiC transport chan- 
nel. We find that heat dissipation is suppressed at room temperatures by the ineffective 
phonon transmission 'T' (lower left panel) across the SiC/metal or SiC/Si interfaces 
and by the phonon Knudsen effect 'K' (lower right panel), that is, the suppression of 
the thermal conductivity within the SiC by interface scattering at layer thicknesses 
c^iayer < 0.5 fim. The phonon transport is calculated analytically at general probabil- 
ities, p+{-) and 1 — p+{-), for specular and diffusive scattering at the upper (lower) 
interface. 
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Figure 2: Schematics of the effective, atomistic representation in a simple-cubic 
(SC) model [2j adapted to estimate the finite-temperature phonon transport across 
SiC/metal (and SiC/Si) interface (vertical dashed line). The SC volume al^ is chosen 
to contain exactly one physical atom but having the average atomic mass Msic = 
(Msi -|- Mc)/2. A similar approach describes the neighboring media 'me' (metal or Si 
substrate). Force constants, -Fp^. sic(mo) fitted to the longitudinal and transverse sound 
velocities [21 El, hnk the SC volumes within the SiC (metal) while elastic coupling 




connect across the interface jSj. 
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Figure 3: Dramatic variation of the room-temperature thermal conductance ax across 
a SiC/general- media interface specified by key differences in materials properties. A 
strong (weaker) dependence arises with the variation in the ratio of material sound 
velocities c (acoustic impedance Z) which specifies the amount of total internal phonon 
reflection (long- wavelength tunneling probability). 
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Figure 4: Knudsen suppression of the effective in- plane SiC transport Kes by interface 
scattering in the typical layered structures (inserts). The upper (lower) panel reports 
the variation with the thickness rfiaycr of SiC conductance channel at room (liquid- 
nitrogen) temperatures assuming either diffusive (solid curves) or specular (dashed 
curves) scattering at the SiC/Si-substrate interface. The traditional approach based on 
the Fourier law of heat flow (dotted curves) is documented to be inadequate due to the 
micron-scale phonon mean-free path. 



15 



